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Background/Philosophy 

- Elevated Material Behavior 

- Impact on Analysis 

- Multiscale Framework/ V i si on 

Recent Advances 

- Theoretical Modeling/Testing 

- Numerical Integration 

- Material Characterization 



ypicalHighTempera ture Applications 
Demand High Performance Materials 


HTHL/SSTO 


• Complex Thermomechanical Loading \ 

• Complex Material response requires Time-Dependent/Hereditary Models: 
Viscoelastic/Viscoplastic 

• Comprehensive Characterization (Tensile, Creep, Relaxation) for a variety 
of material systems 
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Concorde Orbiter Gen 3 A/C 





Important Phenornenological Observations of Behavior of 
Metals at High Homologous Temperatures (T/T m >0.3) 



Strain-rale sensitivity 

lj* 


,-bi, 


Creep Relaxation 




^ Small strain recovery 
Thermal Recovery Dynamic Recovery 


Small reversed 
inelastic strain 
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® Indications of 

Creep hereditary behavior 

Creep-Plasticity Interactions 
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Classic Reason for Introducing 
Unified Viscoplastic Models 
(e.g., GVIPS Class) 





■Pportantphendmenologican Observations of Behavior of 
HpMetafs at High Homologous temperatures (T/T >0.3) 

Cyclic Behavior 

Stress-controlled 


| RauhEtimBfl rai n 


ShaKL'JOWn 


Shakedown 

Behavior 


Strain-controlled 


Ke action o mean stress 


N'onrtf3a\JitLOn or tllcan stnrSi 






Material Behavior Can Significantly Impact Structural 
Response (e.g. Recovery Mechanisms) 

Arnold et al., “ Creep Buckling of a HttrOK 

Cylindrical Shell Under Variable 
Loading”, Jnl of Eng Mech., ASCE, Vol. 

1 1 5, No. 5, pp. 1 054-1 074, 1 989. 


Applied Compressive Stress/Euler Stress = 0.095 
Normalized Initial imperfection - 0.01 



Normalized radial displacement versus normalized time for variable loading 

histories given in inserts 



j- Segment 


Cooling 

channel 


Unified ViscoptasticModels Capture Deformation 
Response in Rocket Engine Nozzle Liners 

Experiment (GRC) 

SSME Nozzle Liner Geometry £^“ lanL - - 


Actual 


Prediction 

Classical Unified 
(Lockheed) (GRC) 


Severe thermomechanical loading conditions result in 
irreversible strains 


Unified viscoplastic models successfully predict the 
experimentally observed deformation trends 

> Arya and Arnold, AIAA, Vol 30, No. 3, 1992 
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Characterization/V alidation 


Structural Analysis 


Detection Techniques 
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CONSTITUTIVE MODELING 

Structural Mechanics Problem ^ T , , ,, 

Need to concurrently address 

three important and related areas 



Knowledge of the material’s life and constitutive 
behavior is a prerequisite for assessment of 
component performance/reliability 


i) mathematical formulation for the 
accurate multiaxial representation 
GVIPS Classes 

ii) algorithmic developments for the 
updating (integrating) of external and 
internal state variables - PEA User 
Definable Subroutines 

iii) parameter estimation - COMPARE 


This approach allows one to overcome the two major obstacles for practical I 
utilization of sophisticated time-dependent (hereditary) models: 

1) lack of efficient and robust integration algorithms - FEA Linkage issues 

2) difficulties associated with characterization of large number of material 
parameters and appropriate experimental “data content” - COMPARE & 
sensitivities 

iASA Glenn Research Center 



www.mdmc.net 


Traceable in-house relational materials database 


MATERIAL 

^PEDIGREE^ 

Batch number 
Compos! lion 


Theory 


COMPARE 


CALCULATE 

key parameters, 

ANALYZE 

mutllple files. 

STORE 

x original da l a 


STATISTICAL SUMNIAfllES 

OF MULTIPLE TESTS 


ABAQUS 


Large Scale Implementation 

• Integration scheme 

• Multimechanism Constitutive 
Relation 


Mathematical Characterization 

N^f 

M aterialBehavior 


Automatically 
input informal 


Source Code 


Ob ject Code 


GVIPS 


FEA Analysis of 
component 
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The Desired vision For Design and Analysis 









Damagi 


fhermomeChanical Testing in Support of 
I Constitutive Model Development 


Characterization 
Tests 


Exploratory 

Tests 


Identify Fundamental Def & Damage 
Mechanisms 


Illuminate Salient Material Response Features 

• Isotropic/kinematic Hardening 

• Time Dependent/ Time-Independent 

• Sensitivity Hydrostatic Stress Field 

• Isotropic/Anisotropy Material Symmetry 

Guide Mathematical Structure of Model 

Guides Specimen design/ Test Method 
Development 


Validation 
Tests 


Often structural in nature 

Provide prototypical response 
data which is to be compared 
with model predictions 

Ideally provide feedback for 
subsequent model refinement 


SA 
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Experimental Observations 


•Reversibility 

> rate-dependent instantaneous stiffness 

> transient creep/relaxation 

> 

•Theoretical demarcation (Exp. Verified) 
•Irreversibility 

> strain- stress dependent 

> nonlinearity 

> strain rate dependence 

> creep with steady- state 

> relaxation with finite residual state 

> creep/plasticity interaction 

> thermal recovery 

> nonlinear kinematic/isotropic hardening 

• Anelastic recovery during reversal in both 
quasilinear and fully developed inelastic 
regions 







and Irre versibleThreshold Surface 



/~\Time independent ^Time dependent 

^reversible region ^reversible region lH ' |rrevers| bie 

(quasilinear) 
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Experimentally verified for both 
TIMETAL 21S and Ti-6-4 \ I 

GR Cop-84 doesn ’t appear to 
exhibit strong viscoelastic response 





Theoretical/Computational Motivation 

In view of four + decades of active research in the area of inelastic behavior 
modeling, the need still exists for an: 


Accurate representation of material response details over an 
extensive domain of time, stress, temperature, loading conditions ... 


Assessment 


Technical 


Practical Implication 


Non-associative 

- Nonsymmetric Tangent Stiffness 

- Coupled system of Stiff Diff. Eq. 

=> Non-uniqueness of solution 
=> Implementation into large scale FEA codes 
problematic 
=> Difficult to integrate 

Numerous nonphysical material parameters 

=> Requires expertise to characterize model 

Single-mechanism models 

=> Qualitatively capable, yet quantitatively 
limited in response spectrum 
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Utilize ConcepiofTJhermodynamic Internal State 
Variables to Obtain Constitutive Equations 
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Dissipation Potentials 
Complementary Type 
Q = Q (variables) 


Equations of State 

irmodynamic Potentials 


Thermodynamic 

Conjugates 

Displacement-Like 





Dissipation Potentials 
Free Energy Type 
'P = 'F (coniusates) 



Missing Link in past potential 
based theories 


Glenn Research Center 



Provides a consistent framework for deformation and 
damage modeling 

Nonisothermal and/or anisotropic extension straight forward 
Nonproportional loading histories automatic 

Automatic satisfaction of the Dissipation Inequity of 
Thermodynamic s 

Eliminates the “ad-hoc” nature of model development 

Provides sufficiently general variational structure. 

Constitutes cornerstone of regularity and bounding (or 
limit) theorems in plasticity and viscoplasticity. 

Lends itself to robust numerical implementation 
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Physical Mechanisms Underlying The 
Partitioning of Energy ^ Complementary Type 



Equations of State 

◄ N 


<*> = + ^IR 


(e,e‘); Y; (e ve ,p) 

Displacement-Like” 


sXotal = Stored + Dissipated 
oe = O + D 


x 



n = n R 

^IR | 






Stored (<I>) = Reversible + Irreversible 

A 

ice Di 


Lattice Distortion Dislocation Pile-up 
^ ^ 


Reflects change in microstructure 


Dissipation (£2) = Reversible + 

T >efo m i'ii 

Dislocation bowing 






Reflects mobility/rate of evolution in 
microstructure 


Irreversible = (deformation) + 

Glide/plastic Slip 


(diffusional; mass/vacancy) 


Thermal recovery 
Dislocation/boundary interaction 
Formation of cell structure 
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General Multimechanism Hereditary 
Behavior Model of the GVIPS Class 


0_o r +cd ir Q-Q r +q, r 




Specific Choice of Energy Potentials and Material 

Functional Forms 


&R = qjj' 1 q|£ J + ^ Oy ? p{j 
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Dissipation 


where: 
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(C (B ) = %>{l - s/cw} L 


(*h 




Results Illustrating 
Recent Improvements 
Made to the Hardening 
Functional Form in GVIPS 

Model 


Previous Non-saturating 

g(G)=H / G B 

Current Saturating Form 

g(G)=H(1-G) B 

G = [V^CCy OCjj)/ K 2 (b) ]°' 5 


TIMET AL 21 S: 650°C 
Strain Controlled Tensile 
Single Mechanism 
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Stress, MPa Stress, MPa 


500 - 







Comparisdn^fSpecific Hardening 
■ Forms UnderGycilc Loading 




TIMETAL 21S: 650°C Strain Controlled 
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New Saturatiftg-Form Does Not Adversely Impact 
Ability to Represent Creep/Relaxation 




But need at least two 


mechanisms to capture both 
creep and relaxation well 


400 - 


300 - 


of 200 

V) 

£ 

35 

100 H 


0 J 



(c) 


20 40 60 

Time, 1000 sec 


80 
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Robustlnteg ration Scheme Key For 
C icient InelasticFi 


Common approaches for integration of Advantages of Implementation 


rate equations: 

1) Non-Iterative : explicit; semi-implicit 

No local iterations 0 less overhead 
stability problems 

2) Iterative : fully-implicit 


Requir 


overhea 


Unconditional stability 

Consistent Tangent Stiffness 0 

Quadratic Convergence of global 
Newton-Raphson Iterations 


Directly applicable for 3-D and sub- space 
loading(planeStrain, axisymmetric, etc) 

Generalized Material Symmetry Operators 
(which influence flow, hardening, 
recovery, relaxation spbctrum, etc.) 

Efficiency (through explicit algorithmic 
tangent stiffness) 

Robustness (through “slack” Mne search) 


Selected: 

Backward Euler with Line Search 
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Results lllustrating the Efficiency of The Numerical 
■■■■! Implementation of GVPS 


Backward Euler wjth Line Search 


W/Kanthal , f=2xl0 _3 /s , £ =0.0144 


method 

number of 
load steps 

CPU time 

GIT 

LIT 

explicit 

10,000 

180.0 

3 ' 

0 

implicit 

100 

5.0 

2 

4 

implicit 

10 

1.05 

4 

10 

implicit 

5 

1 ( 54 s) 

10 

20 



Strain (%) 


Under cyclic conditions 

W<'' 
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</> 


-150 


^Explicit Failed 


1 00 Steps/path 
• 1 0 Steps/path 

8 

1 

_ \ 0 
\ 1 * 

3VA 

T *■ ^*5 

1 

— • J 

7 

l 1 


0.5 1.0 

Axial stress MPa 


Under nonproportional loading 
conditions 



Involves Sufficient “Data Content 


w 


Viscoelastic Material Parameters 

2+2M number, i.e., E s , v, (M (a) , p (a) ) 

Viscoplastic Material Parameters 

- Flow K, \i, n 

- Hardening H b , and p, 

- Recovery: R b and m b 

3 + 5N irreversible material constants 

Types of Experimental Tests 

Strain controlled Tensile Tests (multiple rates) 
Creep Test (Monotonic and/or step) 

Relaxation (Monotonic and/or step) 

Cyclic Tests (Fully reversed, ratcheting) 
Biaxial Tests (tensile, creep, relaxation, cyclic) 
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Temperature 








Desire a mixture (rather than 
numerous of one type) of teStf at 

numerous temperatures 

■ 

i 

Quality vs. QuantitB 



PARE CORE 
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Sensitivity 

v y 


Direct Differentiation Approach 
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1 

L * 

Analyze 

r ] 
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Sequential Quadratic 
Programming (SQP) 


COMPARE 
(driver) 



mplicit Integration 
for Primal Analysis 


Identify active/passive variables for a test 

Scaling design variables and objective function 

Formulating a single design optimization problem 
weighted objective function. 

Constraints 

sensitivities 





[ F 

Result! 

5 

V 




• Final Optimum Material Parameters 

• Combined & Individual Error Functions 
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Inelastic strain (%) 


Comprehensive Characterization of 
The Deformation Response of TIMET AL21S 



* Exp. 103 MPa 


• Exp. 230 MPa 
53 Exp. 345 MPa 


Pred. 103 MP 3 ; 

Pred. 238 MPa 

- 

Pred. 345 MPa 

In 


^ t 


Experiment 

Prediction 


110 MPa 


- 72.4 MPa 


10 20 30 

Time (1000 sec) 


Loading rate 
68.95 MRa/sEC / A 


a. 200 
2 




tb) 



30 eo 

Time (1000 sac) 


^ o 


0.2 0.4 0.6 

Total strain (%) 


^400 




Total strain rate sec 
5x1 CH 


Step 3 


-1.50 -0-75 Q.OT 0.75 

Total fiSfatn [%} 



Wide Range of Application 


Stress: l->60Ksi 
Time: 2 ->90000 sec 
Temp: 650 C 

Loading Rates: 10" 2 ->10" 10 


i 



Characterizatiojiof IN738LC @ 850 # C 

Loading 

Elastic + 4 Viscoplastic Mechanisms 






Final characterized parameters using four Viscoplastic 
mechanisms for IN738LC @850° C 

Material 

Parameter 

Units 

Value 

Material 

Parameter 

Units 

Value 

E 

MPa 

1.5xl0 5 

Pi 

- 

1 ( 6 )* 

V 

- 

0.33 

^2 

- 

1 ( 6 )* 

K 

MPa 

0.1 

P3 

- 

1 ( 6 )* 

Ki 

MPa 

61.43 

P4 

- 

1 ( 6 )* 

k 2 

MPa 

64.37 

Ri 

1 /s 

1.0x1 O ' 21 

k 3 

MPa 

62.30 

r 2 

1 /s 

1.0x1 O ' 21 

K 4 

MPa 

75.08 

r 3 

1 /s 

l.OxlO 21 

n 

- 

1.486 

R4 

1 /s 

l.OxlO 21 

ft 

MPa -s 

3.79xl0 14 

Hi 

MPa 

4.6xl0 4 

mi 

- 

0.001 

h 2 

MPa 

5.13xl0 4 

m 2 

- 

0.001 

h 3 

MPa 

8.33xl0 v 

m 3 

- 

0.001 

H 4 

MPa 

9.458xl0 7 

m4 

- 

0.001 





* the value between parentheses was determined in the FE simulation of the 
experiment 
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% £/W 




a 22 (MPa) 


Multiaxial stress relaxation 



0 20 40 60 80 100 120 140 160 180 

time (sec) 


Multiaxial ratchetting 



Number of Cycles 


Correlation of GRCop-84 Utilizing 
Multimechanism GVIPS Model 


1 VE mechanisms 
4 VP mechanisms 


Relaxation Tests 




4* 


A ^ ^ < 


20000 40000 60000 

time, sec 


80000 100000 




Creep Tests 
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Structural Verification Testing 


• Ideally should provide 
feedback for subsequent 
model refinement 

• Provide prototypical 

response data which is= 
to be with 

model predictions 





Consequently: ^ m 

• Need accurate temperature, strain and load information at a 
variety of locations - required for any true validation I 

• Number of cycles to failure (alone) not enough 

• Instrumentation incredibly challenging (sever environment) 
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SummaiydfiAdvances in Material Modeling 
H ( SynergisticTechnology) 


Generalized, Fully Associative, 
Multimechanism, Viscoelastoplastic 
Model Available 

Reversible/Irreversible Regimes 

- Spanning wide time, stress, 
temperature spectrum 

- Nonlinear Hardening with 
Saturation 

Ability to capture ratcheting 
Stiffness and/or Strength Reduction 

Automated Material Model 
Characterization 

- via COMPARE 
Materials thus far: 

• Ni based; Cu based; Ti 

• MMC and PMC 
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licit Integration Algorithms 

Directly applicable for 3D/sub-space 
loadii 



Generalized Material Symmetry 
Operators (which influence flow, 
hardening, recovery, relaxation 
spectrum, etc.) 

Efficiency (through explicit 
algorithmic tangent stiffness) 

Robustness (through “slack” line 
search) 


\ 


• Now Commercially Availabk 
- COMPARE 

GVIPS - via UMATs 



OoerTCh 


pen Channel Software 


www.openchanneIfeundation.org 



Experiments 




37446 . 27292 
37448 .29514 

37450.27153 
37453.27222 
37455.28194 
37457 .26597 
37460.45486 
37461 . 67014 
37462 .34931 
37464 .27778 
37467 .28403 

37469.27153 

Data 
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0.838069 

1.320609 

1 .793447 

1.005082 

1.068973 

1.040952 

1.268637 

1.357503 

1.089025 

1.13265 

1.096359 

1.064865 


COMPARE fits the 
GVIPS material 
parameters to 
experimental data 
within minutes. 


Q 

t 

© 

t 

\ 

1 

© 

GVIPS 


User 

Definable 

Material 

Model 


The resulting UMAT can 
be immediately 
accessed by the Finite 
Element Analysis 



Finite 
Element 
Analysis 




■uture Work 


Extend formulation to account for 

> Coupled Nonisothermal Issues 

> Probabilistic Material Behavior 

Characterize additional material systems 
Verify under prototypical loading histories 


Implement softening (damage) mechanisms into 
COMPARE - theory complete 

> Characterize strength/stiffness reduction parameters to 
account for softening effects 
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Questions? 
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